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Abstract 

We propose a three loop radiative neutrino mass scenario with an isolated doubly charged singlet 
scalar AA 1 * 1 without couplings to the charged leptons, while two other singly charged scalars hf 
and hf attach to them. In this setup, the lepton flavor violation originating from AA 1 * 1 exchanges 
is suppressed and the model is less constrained, where some couplings can take sizable values. 
As reported in our previous work |jlj], the loop suppression factor at the three loop level would 
be too strong and realized neutrino masses in a three loop scenario could be smaller than the 
observed minuscule values. The sizable couplings can help us to enhance neutrino masses without 
drastically large scalar trilinear couplings appearing in a neutrino mass matrix, which tends to 
drive the vacuum stability to become jeopardized at the one loop level. Now the doubly charged 
scalar AA^ has less constraint via lepton flavor violation and the vacuum can be quite stable, and 
thus a few hundred GeV mass in AA^ is possible, which is within the LHC reach and this model 
can be tested in the near future. Note that the other hf and hf should be heavy at least around a 
few TeV. We suitably arrange the charges of an additional global JJ{ 1) symmetry, where the decay 
constant of the associated Nambu-Goldstone boson can be around a TeV scale consistently. Also, 
this model is indirectly limited through a global analysis on results of the LHC Higgs search and 
issues on a dark matter candidate, the lightest Majorana neutrino. After hf and hf are decoupled, 
this particle couples to the standard model particles only through two charge parity even scalars 
in theory and thus information on this scalar sector is important. Consistent solutions are found, 
but a part of them is now on the edge. 
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I. INTRODUCTION 


Recently, the second round of the physics run at the Large Hadron Collider (LHC) started 
and a magnificent operation for exploring the scalar sector describing the electroweak (EW) 
scale was launched. The greatest achievement at the first run of the LHC is the observation of 
the new Higgs-like scalar boson, which was the last missing piece of the standard model (SM), 
around 125 GeV by the ATLAS and the CMS experiments 2, 3J]. 

The Yukawa couplings to heavy fermions of the SM, namely the top and bottom quarks 
and the tau lepton have been surveyed with good precision, whereas the lighter states are 
still mysterious from both the experimental and theoretical point of view. The extremely 
tiny observed masses of active components of the neutrinos would be key to investigating 
the scalar sector theoretically because we should accept at least 10 12 -order hierarchy in the 
Higgs Yukawa couplings for describing the neutrino nature within the SM. 

One of the most stimulating ideas for addressing this issue is radiative generations of the 
neutrino profiles. Loop suppression factors should appear in the neutrino masses in this type 
of scenario, which help to alleviate the hierarchy in couplings. Another motivation for this 
direction is that the continuous and/or discrete symmetries ensuring the loop origin could 
also guarantee the existence of a (or multiple) dark matter (DM) candidate(s). Following the 
andmarks jfi-sj], recently a variety of works on radiative seesaw model have been done jfl- 


591 ]. where we can also find studies emphasizing non-Abelian discrete symmetry 60-6J , 


radiative generation of quark/charged lepton masses [68-72], operator analysis [73|, 74, 
radiative models accompanying conformal EW symmetry breaking 75H77], and others 178 - 
90J. 


From the naturalness point of view, higher-loop generation is better. The first three loop 


model for a natural explanation of the neu 
following works continue jl, 
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es was proposed in Ref. ]7] and the 


, where in such situations, 


couplings related with the neutrino masses can be close to unity compared with those in 
models with one or two loop level generation. 

On the other hand, three loop generations could face problems owing to the largeness 
of couplings. As discussed in Ref. jlj, the three loop suppression factor 1/(47t) 6 sometimes 
looks very strong and we probably enhance part of couplings for a suitable realization of 
the neutrino masses. Considering large Yukawa-type couplings with lepton flavor viola- 
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tion (LFV) tend to result in unacceptable enhancements in LFV processes. Hence, choosing 
sizable scalar trilinear couplings appearing in the neutrino mass matrix, which do not gen¬ 
erate LFV directory, seems to be a reasonable prescription. But this diagnosis could still 
be at least partly a misjudgment. When we go for one loop level, as shown in Ref. jl|, 
these substantial trilinear couplings give drastic negative contributions to quartic couplings 
of charged scalar bosons and the vacuum can be threatened with being destabilized. The 
last option would be to accept decouplings of additional charged bosons, namely around 
10 TeV, where the neutrino profile itself can be suitably generated. But finding clear signals 
at colliders becomes very difficult even at the updated 13 or 14 TeV LHC. 

To circumvent the situation, we propose a refined three loop model by use of a global 
17(1) symmetry without additional discrete symmetry. A key point is introducing additional 
Majorana neutrinos, which can violate lepton flavors in the fermion line inside diagrams 
describing the neutrino masses. It is important that with the above setup, a doubly charged 
scalar fc ±=b no longer needs to have direct couplings to the charged leptons. We mention that 
in our previous model without Majorana neutrinos in Ref. [1], like in the Zee-Babu model js], 
a doubly charged scalar should attach the charged leptons to generate violation in lepton 
flavors, where tree-level lepton flavor violating processes are generated and we cannot put 
large values in the corresponding couplings consistently. 

In the present k^-isolated scenario, the doubly charged scalar k ±zL is quarantined from 
the charged lepton sector at the leading order by a suitable choice of global 1/(1) charges. 
Now, /c ±=b cannot contribute to phenomena with LFV at the leading order and constraints on 
couplings are weakened. Consequently, the scalar trilinear couplings can take smaller values 
and the vacuum stability would not be so serious even when we consider a few hundred GeV 
fc ±=l= , the detection of which can be an evident signal for probing this scenario at the LHC 
experiments. 

This paper is organized as follows. In Sec. [Hi first we introduce our setups and subse¬ 
quently, we discuss miscellaneous issues in this model, namely, forms of scalar masses and 
mixings, properties of the Nambu-Goldstone (NG) boson associated with the breakdown of 
the global 17(1) symmetry, sizable correction via charged scalars to vacuum stability, form 
of the active neutrino mass matrix at the three loop level, and details on processes accom¬ 
panying LFV in order. In the following Sec. IIIII after having a discussion on analogies with 
the Zee-Babu model jf]], we execute parameter scans both in the normal and inverted hier- 
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archies. In Sec. 1333 we make global fits of signal strengths of the Higgs production in various 
channels announced by the ATLAS and the CMS experiments, which restrict possible values 
of the mass of the doubly charged scalar and the mixing angle between the SM Higgs boson 
and an additional charge parity (CP) even scalars, whose vacuum expectation value (VEV) 
breaks the global 17(1) symmetry. In Sec. [Vj we discuss properties of the dark matter can¬ 
didate of this scenario, which is the lightest right-hand neutrino, through a relic density 
calculation, an estimation of the invisible decay width of the observed 125 GeV scalar and 
an evaluation of spin-independent direct detection cross section. Section [VT] is devoted to 
summarizing results and making conclusions. In Appendix [A] we give analytic forms of the 
loop functions describing lepton-flavor-violating processes. In Appendix |B] a part of partial 
decay widths of the two CP even scalars with nontrivial forms is described. In Appendix Id 
we summarize the averaged matrix elements squared for relic density calculation. 


II. BASIC ISSUES ON THE SCENARIO 

A. Model setup 

We discuss a three loop induced radiative neutrino model with a 17(1) global symmetry. 
The particle contents and their charges are shown in Table [B We introduce three Majorana 
fermions Nr 123 and new bosons; one gauge-singlet neutral boson E 0 , two singly charged 
singlet scalars (hf,hf), and one gauge-singlet doubly charged boson fc ±=t to the SM. We 
assume that only the SM-like Higgs $ and the additional neutral scalar E 0 have VEVs, 
which are symbolized as ($) = v/\f2 and (E 0 ) = v'/\/2, respectively. Here, we set v as 
~ 246 GeV. It is natural that v' is greater than v to some extent to evade a large mixing, 
which is constrained by the LHC data; see Sec. HV] for details, x (^ 0) is an arbitrary 
number of the charge of the hidden U(l) symmetry, 1 and under the assignments, neutrino 
mass matrix is generated at the three loop level. Note that a remnant Z 2 symmetry remains 
after the hidden 1/(1) symmetry breaking and the particles Nr 123 and hf have negative 
parities. Then, when a Majorana neutrino is the lightest among them, the stability is 
accidentally ensured. 

1 Here, we assume that this 17(1) symmetry is global. However, one can straightforwardly move to the local 
one by introducing additional fermions 
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Lepton Fields 

Scalar Fields 

Characters 

L U 

e Ri 

Nr, 

d> 

So 

ht 

ht 

k ++ 

SU(3) C 

1 

1 

1 

1 

1 

1 

1 

1 

SU(2) l 

2 

1 

1 

2 

1 

1 

1 

1 

U(1) y 

-1/2 

-1 

0 

1/2 

0 

1 

1 

2 

U( 1) 

0 

0 

—x 

0 

2x 

0 

X 

2x 


TABLE I: Contents of lepton and scalar fields and their charge assignment under SU(3)c x 
SU(2)l x U(l)y x U{ 1), where 17(1) is an additional global symmetry and x / 0. The sub¬ 
scripts found in the lepton fields i (= 1,2,3) indicate generations of the fields. The bold letters 
emphasize that these numbers correspond to representations of the Lie groups of the NonAbelian 
gauge interactions. 




global symmetry is given 


as 


-C 


Y 


1 1 

= {ye)ijLL^eRj + -( yL)ijL c Li L Lj hi + (yn)ijN Ri e c Rj h 2 + -(y N )i j T, 0 N^ t .N Rj + h.c. 


V = m ||$| 2 + m||E 0 | 2 + m 2 hl \ht \ 2 + m 2 h 2 \hi \ 2 + m 2 k \k 
+ XuT,lh^h^k ++ + ^ 22 ^ 2^2 ^ + h.c. +A$|$| 4 - 

+ A$h 2 | < h| 2 |/b /| 2 + 


(III) 


?IL ++ I 2 

A$|$| 4 + A$ s |<f>| 2 |E 0 | 2 + A^J'h \hl\' 


AnS 0 ft 1 

+ A$fc|$| 2 |/c + 2 + As|E 0 | 4 + As/n |E 0 | 2 |/i^| 2 + A y/, 2 
, A Sfc |E 0 | 2 |fc ++ | 2 + A hl |h+| 4 +A hl/! . 2 |h+| 2 |/i+| 2 + X hlk \ht\ 2 \k 
+ Xh 2 \htY + x h 2 k \h 2 \ 2 \k ++ \ 2 + X k \k ++ \\ 


|S 0 |>+| 2 


2 lh++l 2 


(11*2) 


where the indices i,j indicate matter generations and the superscript “c” means charge 
conjugation. 2 y R , y R and y n are antisymmetric, general, symmetric three-by-three matrices, 
respectively. The first term of Cy generates the charged-lepton masses following the SM 
manner. Majorana mass terms are derived from the fourth one after So obtains a VEV. 
Note that this VEV also generates an effective trilinear coupling fin (in front of h/h//c ++ ) 

2 For SU(2)l doublets, charge conjugation is defined with the SU{2)l rotation described by a Pauli matrix 
as to 2 . 
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from the quartic coupling An (in front of £q h 1 h l k ++ ^j , where the coefficient is given as 

An = An(s;) = A„-h, (II.3) 

where we use the parametrization of (So) declared in Eq. (III.5(1 in the next subsection. We 
assume the following two things: (i) An and H 22 are positive real; (ii) y^ is diagonal and 
obeys the hierarchy (|/v)n < ( Vn )22 < ( Vn )33 among positive-real parameters, which means 
that a generated Majorana mass matrix for Nr is also diagonal one and the mass ordering 
is M Nl < M N2 < M Ns . The concrete forms of the masses are 


M Nl — ^=(r/Ar)n, M N2 — —=(y N ) 22 , M Na — —=(y N ) 33 . 


V2 




V2 y 


(n.4) 


B. Mass eigenvalues and eigenstates of scalars 


The neutral scalar fields are parametrized as 


<f> 


v+(f>-\-iz 


£ 0 = 


V _l±^ e iG/v' 

V2 


(n.5) 


where v ~ 246 GeV is the VEV of the Higgs doublet field, and w ± and 2 are (would- 
be) NG bosons that are absorbed as the longitudinal components of the W and Z bosons, 
respectively. Requiring the tadpole conditions, dV/d(p\ v = 0 and dV/da\ v / = 0, the resultant 
mass matrix squared of the CP even components (0, a) is given by 


m 2 (0,a) 


2A$u 2 \$y,vv' 
X^’evv' 2X^v' 2 


cos a sin a 
— sin a cos a 



m 2 h 

1 

O 

1_ 


0 

m H _ 



cos a — sin a 
sin a cos a 


(II.6) 


where h is the SM-like Higgs (rrih = 125 GeV) and H is an additional CP even Higgs mass 
eigenstate. The mixing angle a is determined as 


sin 2a = 


2Xq,Y J vv' 


m 


H 


mi 


(IL7) 


The neutral bosons 0 and a are rewritten in terms of the mass eigenstates h and H as 


0 = h cos a + H sin a, a = —h sin a + H cos a. 


(US) 


A NG boson G emerges due to the spontaneous symmetry breaking of the global U{ 1) 
symmetry. The mass eigenvalues for the singly charged bosons hf, hf and the doubly 
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charged boson /c ±d= are given as 


m 2 h ± = m hi + + X^ hl v' 2 ), m 2 h ± = m 2 h2 + -(A ^> h2 v 2 + A m 2 v' 2 ), 

m 2 k ±± = m\ + i(A 4 > k v 2 + A Sfe n' 2 ), 


(IL9) 


where the three charged particles are not mixed due to the symmetries of the system and 
thus they themselves are mass eigenstates. 

C. Issues on the Goldstone boson 


Accompanying the spontaneous breakdown of a U( 1) global symmetry, a NG boson 
emerges as an almost massless state in theory, which could play significant roles in particle 


9l|. Like the usual Majoron case 92], our NG boson G communicates 


physics and cosmology 
with the lepton sector. 

An important characteristic of G is that, as described in Table |U the lepton doublets 
and the charged lepton singlets do not hold nonzero charges of the global U(l). This means 
that no anomaly-induced interaction to two photons is generated in our setup, which puts a 


significant constraint on the decay constant of NG bosons 


93]. Thus, in the present scenario, 


we can choose “lower” values around a TeV scale without doing any harm. 

Another route for constraining models via NG bosons is through the active-sterile com¬ 
ponent mixing as through the lepton-flavor-violating transition like fi~ —> e~ (NG) seen in 
Majoron seesaw scenarios, e.g., discussed in Ref. [94]. Different from such a situation in our 


case, the active and the sterile components cannot mix with each other since this mixture 
is prohibited by the residual Z 2 symmetry after the global 17(1) breaking shown in Table [U 
Then the absence of this type of constraint is assured via the accidentally remaining sym¬ 
metry. The neutrinoless double beta decay via W exchanges does not restrict our scenario 
since the sequence with W boson also requires the active-sterile mixing. Note that the three 
additional charged singlet scalars have no direct coupling to the quarks and are therefore 
ineffective. 


In contrast, the NG boson G couples to the corresponding CP even component a, which 
should mix with the Higgs component of the doublet $. This means that G can contribute 
to physics associated with the CP even scalars. As we see in Sec. [V] the pair annihilation 









process of the dark matter candidate Nr ± is just an example. 3 


Finally, we briefly comment on possible bounds from cosmological issues. For example, 
an effect on cosmic microwave background via cosmic string generated by the spontaneous 
breakdown of the global Ui 1) symmetry possibly put a constraint on our scenario. The 


bound discussed in Ref. 


96| can be interpreted as v' < 10 15 GeV, and thus this issue is 


harmless. On the other hand, as we discuss later in Sec. 1111B1 at least part of scalar 
self couplings tends to be 0(1) (at around the EW scale) owing to the requirements via 
coexistence of the observed active neutrino profiles and the null observation in lepton-flavor- 
violating currents. This trend would lead to blowups of the self couplings a little bit above 
the lower scale that we focus on in this paper. Then it might not be so fruitful to discuss 
issues originating from physics at a higher scale. 


D. Vacuum stability against charged scalar trilinear couplings 


Vacuum stability has to be especially assured in the Higgs potential against contributions 
from electrically charged bosons (hf, hf, /c ±=t ). However, our model has some loop contribu¬ 
tions to leading-order values of these quartic couplings via the scalar trilinear couplings /in 
and fi 22 - When they are large, we should examine the vacuum stability against the effect. 
Here, we examine this issue at the one loop level. Let us describe these quartic couplings as 


follows, 

0 < K = 4 ? + a £>, ( 11 . 10 ) 

o < V, = (11.11) 

0 < \ k = Af + A), 1 ’, (11.12) 

where the upper indices denote the numbers of the order in loops, and the one loop contri¬ 
butions are given as 

4? = -%n| i Fo(m h ±,m k ±±), (11.13) 

X h 2 = -% 2 2 | 4 F 0 (m^±,m fc ±±), (11.14) 

X k ] = — 4|p,n| 4 F 0 (m ft ±, m^±) — 4\[x 2 2\ 4 F 0 (m h ±, rn h ±), (11.15) 


3 Another interesting topics is collider searches for a NG boson through invisible channels (subsequent 
decays from CP even scalars) 95]. 
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with 


F 0 (m a ,m b ) = 


(4vr) 


dxdy5(x + y — 1)- 


xy 


1 m a lo S 


+ yml ) 2 

log 


- 2 ) + ml 


+ 2 


= < (4tt) 2 

1 1 


(ml — ml) 3 


(for m a ^ m b ), 
(for m a = m 6 ), 


(11.16) 


. (4tt) 2 6m^ 

where the form of /in is shown in Eq. (111.3(1 and each of m± and m 2 in F 0 represents a mass 
of propagating fields in the loops. We include these constraints in the numerical analysis 
later. To avoid the global minimum accompanying charge breaking, the following condition 
should at least be satisfied: 


I/X22 1 < VA [ml + m 2 hl + m 2 h2 + m 2 k + m |] 1/2 , A = ^ A*, 

z=all quartic couplings including An 

(11.17) 

where we assume the simplified configuration, r = |$| = \hf\ = |h^"| = \k ++ \ — |E 0 | and the 
above inequality comes from the requirement that r does not have a finite nonzero value. 
The summation is taken over the coefficients of the 17 quartic terms in Eq. (III.2D including 
Eg h^h^k ++ and its Hermitian conjugate. When all of these quartic couplings are assumed 
to take the same value A, the above condition is rewritten as 

M < 4\/X [m 2 hl + m 2 h2 +m 2 k + m|] 1/2 , (11.18) 

where the contributions via m\ and A$ are insignificant and thus neglected. 


E. Neutrino mass matrix 


A Majorana neutrino mass matrix m v is generated at the three loop level via the diagram 
shown in Fig. [1} which corresponds to the coefficients of the effective operators, —\(m v ) ab x 


( u L a ) c VL b - The form of (m u ) ab is evaluated by a straightforward calculation as 

(mu) ab =j^Y- T7I [(yL)aiT r ne i (yR)ik(M Nk )(yjt)k j m£ j (yfyjb] 

^ ’ i,j,k =1 1V k 


x F, 


' m ht m ht ml mj. m n 2 m 2 k±± ' 
M 2 k ’ M 2 k ’ M 2 ’ M 2 ’ M 2 k ’ M 2 k 


(11.19) 
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FIG. 1: A schematic description for the radiative generation of neutrino masses. 


where the mass scale M k = rnax[m h ±, rri h ±, rri £ t , my, M Nk , rn k ±± is used for factorizing the 
loop function F\ as a dimensionless variable. Here, we take the relationship in masses in 
our setup, my, my < Mjv fe , into consideration and then M k has only the index k. F\ is 
symbolically calculated as follows: 

1 P 


F 1 (X ll X 2l X 3l X 4 ,X 5l X e ) = I dxi- 


( 11 . 20 ) 


(A 2 ) 2 (A 3 ) 2 ’ 

pi pi pi 

dX = / dxdydz8(x + y + z — 1) / dad/3d^d8 8(a + /3 + 7 + 5 — 1) / dpdadu 8(p + a + u — 1), 
Jo Jo Jo 


with 



( 11 . 21 ) 

y(y-l)+z(z-l) + 2yz, 

( 11 . 22 ) 

(aY + 8) 2 — 8 — aY 2 - aX, 

(11.23) 

pA (Xi,X 2 , X 3 , X 5 , X 6 ) - aX 4 - cjX 1} 

(11.24) 


A(X 1 ,X 2 ,X 3 ,X 5 ,X 6 ) = 


a((x + y)X 2 + zX 5 ) 


((aY + 8) 2 — 5 — aY 2 — aX)(y(y — 1) + z(z — 1) + 2 yz) 

(3X 4 + 7 X 3 + 8X 6 


X = 


({aY + S) 2 — 5 - «F 2 - aXy 

y(y- 1) 


y 


Y = 


y 


(11.25) 


(11.26) 


,y + zj y(y - 1 ) + z(z - 1) + 2yz’ y + z 

Here, note that the shape of F\ is completely the same as that in Ref. |l| except for the 
content of X 5 in Eq. (III. 191) . 
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Neutrino mass eigenstates and their mixings are evaluated by reflecting on similarities 


87] . The structure of the fermion line is similar to the Zee-Babu 


to the Zee-Babu model 

model 6 ], that is, a rank two model of the neutrino mass matrix with a massless eigenstate 
due to the antisymmetricity of y l- Let us describe the neutrino mass matrix as 

3 


,)ab = (77pMNS^ laS UpMNs)ab = 


^ 1 {.Ul) ai^ij (jJL ^jbi 
i,j=l 


un = ^2mt i (y'fi L )ik(CkM Nh )(y R ) kj rne j , 


k=1 


_ ynH22 ^ 


m h+ m h+ 

—2 _ -2. 0 0 

Ml ’ Ml ’ ’ 


M Nk 2 m 2 


fc±± 


M 2 ’ M 2 


(11.27) 

(11.28) 
(11.29) 


(4vr ) 6 M 4 

where i, j, k run over 1 to 3, m)] iag = (mi, m 2 , m 3 ) are the neutrino mass eigenvalues, and 


E, 98 


7/pMNS (Pontecorvo-Maki-Nakagawa-Sakata matrix 
nalize the neutrino mass matrix, which is parametrized as 87] 


is the mixing matrix to diago- 


fTpMNS — 


C 13 C 12 


C 13 S 12 


-C23S12 - S 23 S 13 C 12 e lS C23C12 - S23Sl 3 Sl2e*' 5 


si 3 e iS 
S23C13 


S 23 S 12 — C 23 Si 3 Ci 2 e 


iS 


iS 


-S23C12 — C2 3 Si3Si2e C23C13 


1- 

1 

0 

1- 

O 


0 e^ /2 

0 


0 

0 

1 


(11.30) 


with Cjj = cos 9ij and = sin 0 t] . 6 and 0 represent the Dirac CP phase and the Majorana 
one, respectively. Here, we treat the two ratios in F \, m^/My. and m^./Mk as zero values 
approximately, which means that the loop factor has no dependence on i or j as ( 0 . 

Requirements for the observed neutral profiles depend on the ordering of the neutrino 
masses, which are normal [mi (= 0 ) < m 2 < m 3 ] or inverted [m 3 {= 0 ) < mi < m 2]. 4 


When we consider the norma 
observed neutrino profiles 


87j, 


ordering, the following relations should hold for realizing the 


(Vl) 13 = (Sl2C 23 /(Cl2Cl 3 ) + s 13 s 2 3e l5 /c 13 ){y L ) 23 , 
(Vl) 12 = (S12S23/(C12C13) “ Si 3 C 23 e l5 / C\ 3 ) (t/Z,) 2 3, 


{y 


L) 23 


CU33 ~ ?™3 C 13 S 23 + m 2 ^{Ci 2 C 23 ~ ^S12S13S23) 2 , 


iS , 


(y. 


L) 23 


U>23 « -m 3 c 2 13 c 23 s 23 + m 2 e l,p (C 12 S 23 + e ld C 23 Si2Si3)(ci 2 C23 - e ld Si 2 Si 3 S 2 3), 


iS , 


AS , 


4 More details can be given in Ref. 87( for both of the cases. 


12 


























(Vl)23 

2 


(11.31) 


+22 ~ m 3 ci 3 C2 3 + m 2 e l<t) (ci2S23+e ld C23S 1 2S 13 ) 2 , 


id , 


where we use In the case of the inverted neutrino mass hierarchy, the condi¬ 


tions are deformed 


87j, 


(Dl) 23 


(l/L)l3 ~ — (Cl3-S23e tS /Sl 3 ){yL)23, 

(Vl)i2 = +(Cl3C236 l& /S 13 ) (2/l) 2 3, 

2 

W33 ~ ^i(c 23 Si 2 + e* 5 Ci 2 Si 3 S 23) 2 + m2e l<t> (c 12 C23 - e lS S12S13S23) 2 , 


(y 


L) 23 


w23 ~ "ll(Sl2S 2 3 - e l<5 Ci2C23Sl3)(c 23 Sl2 + e l5 Ci 2 Si 3 S 23 ) 

+ m 2 e l<t> (ci2S 2 3 + e*' 5 c 23 Si2Si 3 )(ci2C 23 - e*' 5 Si 2 Si 3 s 23 ), 


( 2 /. 


L)23 


UJ22 ~ mi(si 2 s 2 3 - e +i<5 Ci2C2 3 Sl 3 ) 2 + m 2 e J0 (ci2S23+e*‘ ) C2 3 Sl2Sl3) 2 - 


id, 


(11.32) 


In the numerical analysis in the following Sec. IIII B1 we adopt the latest values of 9 12 , 9 23 , 
013 , Amji_ and A m 3i (l — 1 for the normal hierarchy and i = 2 for the inverted one) found 


in Ref. 


991 ]. The two CP-violating phases <5 and </> are treated as free parameters." 


F. Lepton flavor violations and the universality of charged currents 


In the present model, owing to interactions containing charged scalars, new contributions 
to several lepton-flavor-violating processes are found at the tree or the one loop level. Also, 
universalities of charged currents are threatened by the vertices. Some of them are very 
similar to the case of the Zee-Babu model discussed, e.g., in Ref. |81], while situations are 
changed in the others. Here, we briefly summarize important points: 


• £i —» iji/u: like in the Zee-Babu model, the processes receive additional contributions 
that are never found in the SM via the exchange of the singly charged scalar coupling 
to the charged leptons hf, which have an influence on the value of the Fermi constant 
via charged-lepton decay. If the masses of the right-hand neutrinos are small, there is 

5 The best fit values are as follows: sin 2 6*12 = 0.304, sin 2 6*23 = 0.452, sin 2 0\ 3 = 0.0218, Am^ = 7.50 x 
10~ 5 eV 2 , A m 31 = +2.457 x 10 _ 3 eV 2 (in the normal hierarchy); sin 2 $12 = 0.304, sin 2 023 = 0.5 79, 
sin 2 0i3 = 0.0219, A= 7.50 x 10~ 5 eV 2 , Am| 2 = —2.449 x 10 _3 eV 2 (in the inverted hierarchy) [ 99 ]. 
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a possibility that the channels £i —> £jN\N\ via the exchange of hf are open and they 
fake contributions to £i —> £jvv. In this paper, we do not consider such a light N\. 
Then, we simply ignore this effect in the following analysis part. 


• £y —* £j y: like in the Zee-Babu model, these decay sequences (when m e - > m t -) are 
radiatively generated at the one loop level via the diagrams containing charged scalars 
which the photon couples to. In these processes, both of the interactions, Ly.L^.hi 
and (also with their Hermitian conjugations included) make nonzero con¬ 

tributions to the decay widths. The doubly charged scalar k ±zL has no contribution 
at the one loop level due to the absence of tree-level interactions with the charged 
leptons. 


£j —> £~j'£y£y: being different from the Zee-Babu model, absence of direct interactions 
between the charged leptons and the doubly charged scalar /c ±=t generates the situation 
that these processes are induced at the one loop level. Two types of box diagrams 
contributing to the processes are there: containing (i) Ly^Li.hf interactions; and 
containing (ii) N^e^.h^ ones (where their Hermitian conjugations are also included), 
while two types of photon and Z penguin diagrams are also generated. Here, the 
contribution from the penguins is directory related to the dipole diag ram s of £y —> f?/y 
in the language of branching ratios, e.g ., when we consider /i —> 3e 


100 


102] as 


Br(/i —y 3e) ~ 


«EM 

37T 


log ( —y 

V mi 


11 

T 


Br(/i ey), 


(11.33) 


where Oem is the electromagnetic fine structure constant and the penguin-type con¬ 
tributions to /i —y 3e are suppressed compared with those to /i —> ey. Since the 
experimental bounds on £~ —> f?Jy and £~ —¥ are roughly comparable, we 

simply ignore the contributions from the penguin-type diagrams in the calculation of 
the constraints in this category. 


Now we show the explicit forms of constraints on the present model. We recast the 
results of the analysis on the Zee-Babu model in Ref. 87] for our case, where the three 


£ — > £y processes (/i _ — > e _ y, t~ — > e _ y, t~ — > n~ y), the four types of gauge universalities 
(lepton/hadron, /i/e, r//i, r/e), and the seven £ —> 3£ processes (/i“ —$■ e + e~e~, r~ —$■ 
e + e~e~, t~ e + e”/i _ , t~ — > e + /i~/i _ , t~ — > /i + e _ e _ , t~ — > /i + e _ /i _ , t~ — > /i + ji~n~ ) are 
included. 
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Process 

(*>/) 

Experimental bounds (90% C.L.) 

Cif 

p _ —>• e _ 7 

( 2 , 1 ) 

Br(/i —> e 7 ) < 5.7 x 10 ~ 13 

1.6 x 10“ 6 

t~ —>• e _ 7 

(3,1) 

Br(r —> e'y) < 3.3 x 10 -8 

0.52 

T~ -> /i “7 

(3,2) 

Br(r —> py) < 4.4 x 10 -8 

0.7 


TABLE II: 
analysis in 


Summary of the coefficient Cif in t —>• £7 processes and experimental data used in the 

871- 


Type of universality 

Experimental bounds (90% C.L.) 

Lepton/hadron universality 

E q =d,s,b\V% P \= 0-9999 ±0.0006 

p/e universality 

G e £ p /Gf p = 1.0010 ± 0.0009 

r/p universality 

Gf p /G^ xp = 0.9998 ± 0.0013 

r/e universality 

Gr xp /GT P = 1.0034 ± 0.0015 


TABLE III: Summary of the experimental data on universality of charged currents used in the 
analysis in [87]. 


• —>■ 7 processes: in this case, the result of recasting is 


V 3 

Z-/a=l 

(y l) af (.V l) ia 

(• h,ah,a + I\ ,a) + ELl 

(dll) af (y n)ia 

2 

16mf ± 

h i 

ELl {h,ah,a + If,a) 



r Cif 
TeV 4 ’ 

(11.34) 


with the loop functions 


f*l — X 


h,a = 


I 2 ,a — 


Ka = 


IL = 


(4tt) 2 , 
1 


dx / dy 


x{2x — 1 ) 


o (x + y)m 2 h± + (1 - x - y)ml a 


i*l—X 


dx / dy 


x(2y - 1 ) 


—X 


(4vr) 2 . 

1 


dx / dy 


x{2x — 1 ) 


o (x + y)m 2 ± + (1 - x - y)M 2 Na ’ 


i*l—X 


dx / dy 


x{2y - 1) 


2,a (4vr ) 2 J 0 Jo (x + y)m 2 ± + (1 - x - y)M 2 ’ 


(47t) 2 36m 
1 5 


2 ’ 
± 


( 47 t ) 2 J 0 J 0 (x + y)m 2 ± + (1 — x — y)m 2 a ( 47 r ) 2 36m : 


2 ’ 
± 


(11.35) 

(11.36) 

(11.37) 

(11.38) 


where we use m Va — 0. Here, we treat the final-state lepton £j as a massless particle. 
Concrete forms of the integrals are summarized in Appendix |A] The dimensionless 
coefficient Cif representing the digits in 87] (before recasting) are summarized in 
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Process 


Experimental bounds (90% C.L.) 

■A-ij kl 

/t - —» e + e _ e~ 

(2,1,1,1) 

Br < 1.0 x 10~ 12 

2.3 x 10” 5 

t~ —>• e + e _ e~ 

(3,1,1,1) 

Br < 2.7 x 10" 8 

0.009 

t~ —>■ e + e _ 

(3,1,1, 2) 

Br < 1.8 x 10" 8 

0.005 

t~ — >• e + fi~ fi~ 

(3,1,2, 2) 

Br < 1.7 x 10" 8 

0.007 

t~ — >• /. L + e~e~ 

(3, 2,1,1) 

Br < 1.5 x 10" 8 

0.007 

t~ —>• n + e~ [i~ 

(3, 2,1,2) 

Br < 2.7 x 10" 8 

0.007 

T~ — > /i + /t _ /i _ 

(3, 2, 2, 2) 

Br < 2.1 x 10" 8 

0.008 


TABLE IV: Summary of the coefficient A^m in l — >• 3£ processes and experimental data used in 


the analysis in 


83]. 


Table HU The factor 16 comes from the difference in the coupling convention of the 
interaction L c L .L Lj hl ( L c L ,L Lj h + in the Zee-Babu model). These decay processes are 
one loop induced ones in both of the models, and thus the loop factor 1/(47 t) 2 in the 
integrals is canceled out in the final form in Eq. (111.34)1 . 


• Gauge coupling universalities: in this category, recasting is just straightforward by the 
replacement —y hf , 


(Vl )12 2 
2 


(Vl) 23 

to 

1 

(vl) 13 

2 

2 


2 


(vl) 13 

2 

(vl) 12 

2 

2 


2 


(Ul) 23 

2 

(z/i)l2 

2 

2 


2 



< 0.007 

< 0.024 


f m h ±\ 2 

\TeV/ 

f m hf \ 2 

\TeVy 


< 0.035 


( m *\ 

G eV) 


< 0.04 


m hf 

TeV 


(lepton/hadron universality), 

(/i/e universality), 

(r//r universality), 

(r/e universality). (11.39) 


The corresponding experimental bounds are summarized in Table IIII1 


£ —» 3£ processes: all of the cases are summarized symbolically as 


| A v + B u \ + \An + Bpj\ — 2 Re [-A^vC^y] — 2 Re [-B/vCw] + — |CW| 


1/2 ^ A 


ijkl 


2 ’ 


TeV 
(11.40) 
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with the effective couplings 

A v = (VLyDikiyLyDji J i,o, 
B u = {yLvDuiyLyDjkJifl, 


An — ^2 (yR)al(yR)bk(yR)ja{yR)ib Jl,ab, 

a, 6=1 
3 

Bn = 'Y2 (yR)al(yR)bk(yR)jb(y\l)ia J\,abi 
a,6=1 
3 

Cn = 2 ^ {yR)al(yR)ak(yR)jb(yR)ibM Na M Nb J 2 ,ab- 


(11.41) 

(11.42) 

(11.43) 

(11.44) 

(11.45) 


a,6=1 


The loop functions J\ A b and J 2)0 ,b are given as 


f*l — X 


J l,a6 
*^2,a6 


(4tt) 2 , 

Wj*. 


dx / dy 


1 — x — y 


xMl + yM% + (1 - x - y)m 2 ± ’ 


f*l — X 


dx dy — 


1 — x — y 


2 > 


(11.46) 

(11.47) 


xM N a + y M N h + (1 - x - y)m h ± 

where the form of J\.q is obtained by the replacements Mjy a —> m Ua (~ 0), MN b —* 
m Vb (~ 0) and m fe ± —)• m h ± as 

111 


J 1,0 ^ 


(11.48) 


2 (47t) 2 m 2 ± 

Also in this calculation, we treat the final-state leptons and £~[ as massless 

particles. Other concrete forms of the integrals are summarized in Appendix [A] Here, 
the dimensionless valuables Aijki indicate the numerical values (in the analysis in 
before recasting) and their concrete values are summarized in Table IIV1 The factor 
1/4 originates from recasting. The loop suppression factor 1/(47 t) 2 should appear in 
Eq. (III.4011 since these processes are generated at the tree level in the original Zee-Babu 
model. 


Finally, we briefly have a discussion on the LFV via k ±zL exchange in our model. When 
we focus on the LFV process accompanying two neutrinos, l~ —y k £ + —y 2(h/)* + £ + —> 
(2£~ + 2u) + £ + , where the first step is one loop induced and the two h/ as intermediate 
states are usually off shell since the mass of h/ is bounded from below, where m h ± should 
be more than a few TeV; see Sec. m for details. The order of the loops is the same as that 
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of £ —> 3L But, the off shellness suppresses the decay width of this LFV process. Therefore 
we can simply ignore this possibility. 


III. VALID PARAMETER REGION WITH NEUTRINO SECTOR 


A. Form of yn 


Different from the Zee-Babu model, y r is not a symmetric matrix in general but 
yft ((Mn) yn (= Yr) is symmetric one. Hence, we can utilize htting relations in the Zee- 
Babu model with modifications. Here, we adopt the following forms in yR [and in yjy as 
shown in Eq. flll.dD ] for simplicity, 


Ur — 


M n = diag ( M Nl , M N2 ,M Ns ), 


(HI-1) 


* * * 

* a b 

* b c 

where a, b , c are arbitrary complex numbers. The correspondence to the factors c Oij in 
Eq. (111.28(1 is as follows: 


«« = mitfYidumi,- 


(IH.2) 


Here, because of the mass hierarchy of the charged leptons <C mn 2 < m^ 3 , the terms 
in ujij including are negligible. Combining this issue with the hierarchy assumed in the 
right-hand neutrinos M^ 1 < Mjy 2 < Mn 3 , we conclude that the elements of yn expressed 
by * in Eq. (IIII.1D do not affect the values of uj t j signihcantly. Therefore in the following 
analysis, we only consider the couplings a, b , c of yn in htting the neutrino profiles. 

Now, we rewrite the relations in Eq. (1111.21) as follows: 

= (Yr) 22 — [a 2 (C2^V 2 ) + fr 2 ((3 Mn 3 )] , 


m h 

<^23 

me 2 mt 3 
^33 


m 


h 


— (^r)23 — [ab (C 2M N2 ) + be (( 3 M N3 )], 

= (Yr )33 — [^ 2 (C2^V 2 ) + C 2 (C3 Mn 3 )] • 


(HI-3) 


As shown in Eqs. flll.31h - (111.32jh when we £x the values of { jjl ) 23 , 5 and 0, the values of 0222 , 
CO 23 , W 33 [and also (|/l)i 2 , {vl) 13 ] are automatically determined through the relations. In each 
scanning in the following section, we pick up a solution on a, b , c of the above simultaneous 
equations. 
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B. Parameter scanning 


Adopting the structure of i/r in the previous subsection, we execute parameter scans to 
search for consistent regions in the parameter space. In this model, k ±zL does not contribute 
to the processes with LFV significantly. Thus, we consider the two possibilities m*.±± = 250 
and 500 GeV, while the other two singly charged scalars hi±,h 2 ± have a few TeV masses. 
Here, we select the four choices: m h ± = m h ± = 3, 3.5, 4, and 4.8 TeV. For brevity, we fix 
the three right-hand neutrinos as follows: = rrih/2, Mn 2 = 5 TeV and Mjv 3 = 10 TeV, 

where we mention that Ni is a “Higgs-portal” dark matter candidate. The mass is assumed 
to be around the 125 GeV Higgs resonant region. Detailed discussion on this topic is found 
in Sec. [V] 

We mention that compliance with the relations in the case of the normal hierarchy in 
Eq. (III.311) or the inverted hierarchy in Eq. (1II.32j) leads to the situation that only the 
parameters /in, P22, 5, </>, (t/l)23 are free to be chosen (in addition to the masses of particles 
that we fixed in the above paragraph). Note that the five matrix elements of i/r, namely 
(2/fl)n, (Ur) 12 , (Ur) 13 , (vr) 21 , and (y R ) 31 , are ineffective in the determination of the active 
neutrino profiles, while their nonzeroness possibly modihes the strengths of the lepton-flavor- 
violating processes significantly. Then we consider nonzero values in the five couplings for 
the sake of completeness. In each scanning, we randomly select values of the ten parameters 
within the corresponding ranges 

P 12 = P 22 (= p) e [ITeV, 20TeV], 6e[0, 27r], 0 g[O,27t], (y L )23 e [-1,1], 

W«6 [-0.01, -0.1] U [0.01,0.1] (for (i,j) = (1,1), (1,2), (1,3), (2,1), (3,1)). (III.4) 

Our definition of an allowed point is a set of parameters where all the following require¬ 
ments are satisfied: 

• Observed values in masses and mixings of the three active neutrinos are generated. As 
we mentioned above, realization of this requirement is equal to the compliance with 
the relations in Eq. (III.311) or (III.321) . 

• Predictions for LFVs do not exceed the bounds shown in Tables UTl (l —> £ 7 ), HU (gauge 
universalities), and|TV](^ —> 3£). Concretely, evading bounds corresponds to meeting 
the inequalities in Eqs. (111.34)1 and (111.39)1 — (111.40)1 . respectively. 
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Tflfe itrb 


F ^2) 

f{ C3) 

Number of allowed points 

500 

GeV 

(4.8 TeV, 4.8 TeV) 

-1.16818 

-11.5428 

2804 

(4.0 TeV, 4.0 TeV) 

-2.14198 

-20.0287 

422 

(3.5 TeV, 3.5 TeV) 

-3.30195 

-29.6525 

89 

(3.0 TeV, 3.0 TeV) 

-5.37518 

-46.0932 

16 

250 

GeV 

(4.8 TeV, 4.8 TeV) 

-1.21466 

-11.9695 

1644 

(4.0 TeV, 4.0 TeV) 

-2.24746 

-20.9444 

190 

(3.5 TeV, 3.5 TeV) 

-3.49141 

-31.2321 

31 

(3.0 TeV, 3.0 TeV) 

-5.74153 

-49.0141 

2 


TABLE V: The numbers of consistent points in the normal hierarchy obtained by parameter 
scanning. F ^ 2 "* and F^ 3 - 1 are the loop functions in £2 and £3 shown in Eq. (111.29(1 . respectively. 
Kindly refer to the main body of this subsection for details of this scanning. 


Fulfilling the requirements on vacuum stability in Eqs. (III. 101) (III.121) and (III.181) . In 
the last condition, we adopt the following modified form with physical masses, 

1/2 


1/^22 | < 4\/A 


™it + + m fc±± 


(III.5) 


where we simply change the parameters m^, rrih 2 , rrik to their physical masses and 
ignore the mass parameter mj. The value of the right-hand side would be near the 
original one and it would be useful for roughly estimating this type of bound. In the 
following analysis, we adopt the initial conditions for the quartic couplings, 


d? = A‘“> = A< 0) = 4n, X = 4n. 


(III.6) 


• Ensuring perturbativity, all the couplings should be equal to or less than An. 

Our result in the normal hierarchy is summarized in Table [V] In each of the eight com¬ 
binations of the charged scalar masses, we randomly scan 10 5 points in the ten parameters, 
where the range we consider is shown in Eq. (IHI.dl) . Apart from the previous work [H, the 
result says that the mass of the doubly charged scalar can be light around a few hundred 
GeV, whereas the other two singly charged ones should be heavy (at least) around a few 
TeV. In this scenario, it is very hard to produce hf and hf in colliders. On the other hand, 
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FIG. 2: Histograms showing distributions of <5, q !>, (m) 23 , |( 2 / 21 ) 22 1 and fj, of the result of a scanning 
with the charged scalar masses m k ±± = 500 GeV, rn h ± = m h ± = 4.8 TeV. The total number of the 
consistent data points is 2804. 


detecting k ±i: could be a reasonable option for probing this model in present and future 
collider experiments. 6 

Next, we look into properties of the parameters in the allowed region. Here, we only show 


6 Lower bounds at 95% CL on m k ±± via 8 TeV data were provided by the ATLAS group in Ref. [103 


as 374 GeV, 402 GeV, 438 GeV when assuming a 100% branching ratio to e ± e ± , e ± ^ =l= , 


pairs, 


respectively. In our model, the decay sequence k ±±L —> i ± is one-loop induced and many parameters 

contribute. In this paper, we skip to have detailed discussions on prospects in colliders. 
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I(j/f)22l 


ICy«)22l 


FIG. 3: Correlations between 0 and 5, (j/l )23 and (( 2 / 77 ) 22 !, and n and | ( 2 / 71 ) 22 1 are shown. The 
data set is the same as in Fig. [2j 


the case of m*.±± = 500 GeV, m h ± = m h ± = 4.8 TeV as an example, where we checked that 
other cases are not so different. In Fig. [21 we display the distributions of 5, 0, (i/l) 23 , | ( 2 / 77 ) 22 ! 
and /i based on 2804 data points of the allowed region as histograms.' The correlations 
between 0 and 5, (i/l) 23 an d 1 ( 2 / 77 ) 22 !, and /i and 1 ( 2 / 77 ) 22 ! are also shown in Fig. [3J In the 
present three loop scenario, the loop factor 1/(47 t) 6 in the neutrino masses tends to suppress 
the realized masses very much. Thus, at least a part of the parameters related with the 
masses would be sufficiently large. In the following, we investigate details. 

• We can see that this model could give some preference to 5 around n. On the other 
hand, 6 not around tt is also realized. Few trends are seen in the distribution of 0. 

• A typical digit of the absolute value of (i/l )23 is in 0.5 ~ 0.6, which would be greater 
compared with other radiative neutrino models in one or two loop level. Since |ct 222 1 — 

7 The value of (ur )22 is complex in general obtained as a part of a solution of the simultaneous equations 
in Eq. (IIIDsl) . 
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IW 23 I — |c ^33 1, the other two nonzero components of i/l, (ul) 12 and (i/l) 13 , tend to have 
almost the same order of magnitude. 

• Like the Zee-Babu model, via the relation jo^l — 1^231 — |^ 33 1, the three components 
of the effective symmetric Yukawa couplings Y# meet the definite hierarchy, 

9 

TTL 711 

|(yk)33| : l(^)2 3 | : \(Y R ) 22 \ -^ ^ : 1, (I1I.7) 

my m r 

which means that the original \(yR) 2 2 \ tends to hold a significantly large value. The 
peak of the distribution is around 9, which is still rather small compared with the 
perturbative upper bound In. Because of this characteristic, the masses of the singly 
charged scalars m h ± and m h ± should be greater than around 3 TeV to circumvent the 
bounds. 

• The common coefficients of the trilinear terms /i should be large at around 14 ~ 
15 TeV to compensate for the suppression factor in the neutrino masses. From the 
perturbativity in the couplings An in the form of the trilinear coupling /z n in Eq. (Ill ,31) 
and Dn found in the masses of the right-hand neutrinos in Eq. (111.41) . v' cannot be so 
small. A reasonable value of v' is 0(1) TeV. 

Finally, we briefly comment on the possibility with the inverted mass hierarchy. In our 
model, like the Zee-Babu model, the components of the active neutrino mass matrix ( m u ) a b 
contain the charged lepton masses, and also the Majorana masses with the assumed mass 
hierarchy < Mn 2 < Mjv 3 . Then the normal hierarchy would be preferable. Within 
our search among 10 5 points in the inverted case, even in the choice of (rrik±±, m h ±, m h ±) = 
(500 GeV, 4.8 TeV, 4.8 TeV), no solution is found. 


IV. CONSTRAINT FROM LHC HIGGS SEARCH 


In this part, we evaluate constraints on the parameter space of scalars in this model by 
use of the latest results of LHC Higgs searches by the ATLAS and CMS experi men t groups. 


First, we describe the method we use for global analysis. Like in the papers 


104 


106], we 


adopt the following form of signal strength of the single Higgs production channel with the 
subsequent Higgs decay to the particles /, 

, i Rr> . r 

(IV.l) 


/V 


^total Br^-g./- 

X J 


<y\ 


SM 

total 


Br 


SM ’ 
h-+f 


23 











Process 

ATLAS 

CMS 

Reference 

h -> 

77 

1 1 7 + 0.27 
i - i '-0.27 

1.12 

± 

0.24 

[107,1121 

h —>• ZZ 

* ->• M 

1 44+ 0 40 

1.00 

± 

0.29 

[108,-112] 

h —»■ WW* 

-»■ Ivlv 

1 D9 +0 ' 23 
i.uy_ 0 21 

0.83 

± 

0.21 

[109, .112] 

h -> 

bb 

n ^+°- 4 

U '°- 0.4 

0.84 

± 

0.44 

[ 110 ,- 112 ] 

h -> 

TT 

I 4 + 0-4 

0.91 

± 

0.28 

[m, 112 ] 


TABLE VI: Summary of the latest LHC Higgs search data as (combined) signal strengths in the 
five Higgs decay patterns (h —> 77 ,h — > ZZ* —>• A£,h —>• WW* —»■ lvlv,h —>• bb,h —>• rf). The 

analyses are based on the data sets accumulated in the LHC first run, whose details of the ATLAS 
and the CMS are 4.5 - 4.7 fb " 1 (7TeV) + 20.3 fb " 1 ( 8 TeV), 5.1 fb ” 1 (7TeV) + 19.7 fb " 1 ( 8 TeV), 
respectively. 


where cr tot ai (<++) represents the total production cross section of a Higgs boson in this 
model (SM), respectively. The Higgs branching ratios to the particles / of the SM and this 
model are discriminated by having or not having the superscript SM. 

Here, note that in our scenario, the observed Higgs boson is a mixture of the SU(2)l 
doublet $ and the singlet So as shown in Eq. (III. 81) and no additional colored particle is 
introduced, which means the absence of new contributions to the Higgs production via the 
gluon fusion process. Since the SM gauge bosons and quarks do not couple to E 0 , the ratio 
of the total cross sections is easily calculated as 


^total 


a 


SM 

total 


cos 2 a. 


(IV. 2) 


Whereas, evaluating the ratio of the branding ratios is rather complicated. First, we look 
at the following decomposition, 


Br h-yf 


Br 


SM 

h-+f 


I\ \ 1 

ifV X rf5 


(IV.3) 


where rj ( SM ^ represents the corresponding Higgs total width. In the Higgs decay, the presence 
of the charged scalars in our model (hf, h f, /c ±=t ) also modifies the partial decay widths of 
the photon-associated decay processes, h —> 77 and h —> Z 7 , in addition to the doublet- 
singlet mixing effect. As we discussed in Sec. IIIIBI the two singly charged scalars tend to be 
very massive at around a few TeV at least. In such a situation, we can completely neglect 
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the effect via the two particles and then only the contribution from /c ±=b is included as a 
loop effect in the following part. 

The ratio of the total widths is rewritten as follows: 


T> SM T/i^SM others . SM T/i—>-77 . t-, SM ^h—tZy . l-V—Mnv 

® r h^-SMothers X -pSM ^ ^^h^yy X pSM ^ 7 X pSM ^ pSM ’ 

^ h—>SM others ^ h— >-77 ^ h—7 ^ h 

(IV.4) 


where T^ inv expresses the Higgs partial decay width to invisible pairs, which is written as 


r h 


V SM 


T/i—nnv r^^GG T ^h^rN^ JV/jj . 


(IV.5) 


Concrete forms of partial decay width are summarized in Appendix [0 


Br 


h^SMothers can be described as 1 — Br ^“ 77 — Br ^“ Z7 very precisely. In the following 


„SM 


SM 


analysis, we use the values Br/A 77 = 2.28 x 10 , Br/7z 7 = 1.54 x 10 3 and P)y = 4.07 MeV 
in rrih = 125 GeV reported by the LHC Higgs Cross Section Working Group 113] - The 


ratios of the partial decay widths are described with the help of the formulas for h —> 77 


and h —> Z 7 in Refs. 


114 


11 


as 


i h^f 

rSM 

1 h-+f 


1 h—>■ 77 
pSM 
/i—>-77 


T/H-Z7 
pSM 
1 h—*Zy 


cos 2 ct (/ e (others in SM)) —> 


b/i—>SM others 
pSM 

others 


cos a + 


cos a — 


1 v 2 _ K*{r k ) _ 

2 m 2 k ±± ^ hkk AV{r w ) + N c Q 2 t AJJ 2 (r t ) 


m 


k ±±L 


(ZQkgzkk) Chkk 


Ap 7 (rfc, A fc ) 


2 


? 


cos 2 a, 


with effective couplings 


Chkk = cosaA$fc — since ( — ) A Sfc , g Zk k — ~Qk 


sw 

c w 


and the form factor 


A Zl - 

^SM - v 


±Q A Z l( \ \ \ AT — ^Qt s w) A Z 1 ( \ 

cot 6wA 1 '{T\y, Aw) + - A,/Jr t , A 

swc w ' 


(IV. 6 ) 

(IV.7) 

(IV. 8 ) 

(IV.9) 

(IV.10) 


Nq (= 3), Qt (= 2/3), Qh (= 2 ), Tf' 1 (= 1 / 2 ), cw and sw are the QCD color factor for quarks, 
the electric charges of the top quark, the doubly charged scalar in unit of the positron’s one, 
the weak isospin of the top quark, and the cosine and the sine of the Weinberg angle 9w, 
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respectively. The loop factors take the following forms, 


= — x 2 [2x 2 + 3x 1 + 3(2x 1 — 1 )f(x x )] , 

AfJ 2 (*^) = 2x2 + ( X_1 - 1 )/(^ _1 )] , 

Al\x) = -x 2 [a; -1 - fix' 1 )] , 

A^ 7 (x, y) = 4(3 — tan 2 9 w )I 2 (x, y ) + [(1 + 2a: -1 ) tan 2 9 W — (5 + 2a; -1 )] h(x, y), 
Af^ 2 ( x i V ) = A (®, 2/) - h[x, y), 

Ap{x,y) = h(x,y ), 


(IV.ll) 
(IV.12) 
(IV.13) 
(IV. 14) 
(IV.15) 
(IV.16) 


with the functions 


h(x, V) = + 77^7 [fix- 1 ) - fiy- 1 )] + [cgix- 1 ) - giy- 1 )] , (IV.17) 


2 ix-y) 2ix-y) 2 


h (x,y) = 


xy 


[fix x ) - fiy *)] • 


(IV.18) 


2 ix - y) 

In the above formulas, forms of the input variables to the loop factors r* and A, are defined 
as fractions by the Higgs boson mass (m^) or the Z boson mass fmz) 


4 TYl 2 4-777 2 

Ti = -j-, Aj = — Y (i = t,W,k). 

V m z 


(IV.19) 

m h m 2 

The two ratios usually take values above one (m^ < 2m;, mz < 2 mf). The two functions 
fiz) and giz) iz = a: -1 or y -1 ) are formulated as 


m = 


g(z) = 


arcsm wz 


log 


l+vV^ 

l—AT 




,-l 


— Z7T 


1 — 1 arcsin yfz 


yl- 




■°g i 


1+Vl—z 1 


-AT 


— 


m 


for z < 1, 
for 2 > 1, 

for z < 1, 
for 2 > 1, 


(IV.20) 


(IV.21) 


where the situation mh < 2m,, < 2m; corresponds to z < 1. 

To estimating the consistent parameter region with the latest results of the Higgs search 
at the LHC, we define the y 2 valuable by use of the signal strength in Eq. (1IV. 1 1) as follows: 


x 2 = 


E 


/=77,ZZ* ,WW*, 
bb,TT (ATLAS) 


h/ /v 


+ 


E 


/=77,ZZ* ,WW *, 
bb,TT (CMS) 


l- l f h-/ 


(IV.22) 


Here, we take the results of the five Higgs decay channels reported by the ATLAS and the 
CMS experiments into consideration, which are h —> 77, h —> ZZ* —> h —y WW* —y 
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FIG. 4: The la (blue) and 2a (red) allowed regions of the global analysis based on the LHC 
data summarized in Table EO From top left to bottom, we choose the parameters (A<j>fc, AsaV) = 
(1.0, v), (7 t,v), (0.1, v), respectively. 


£u£u, h —y bb, h —y TT summarized in Table EH The two hatted symbols fif and <7/ repre¬ 
sent the corresponding central value and error, respectively. We assumed that each of ten 
experimental inputs follows the Gaussian distribution and there are no correlations among 
them. Also, when an error is asymmetric, we adopt its simple average as an input value 
of the corresponding by for analysis. These simplifications are justified for our purpose of 
roughly estimating survived regions of the parameter space of this model. 

In the following analysis, for simplicity, we focus on a reasonable situation in which the 
mass of the DM is around rn,h/2 where the Higgs invisible channel to a pair of the DMs 
is near the threshold and negligible. Detailed discussions on the DM candidate are given 


27 



















m k ±± (GeV) 


m k ±± (GeV) 


FIG. 5: The la (blue) and 2a (red) allowed regions of the global analysis based on the LHC 
data summarized in Table m From left to right, we choose the parameters (A$fc, A^V) = 
(1.0, 5v), (1.0, lOw), respectively. 

in Sec. |V] Now, apparently, we see the five parameters m k ±±, A$fc, Asfc, v' and a govern 
the signal strengths of the Higgs boson. But in fact, as shown in Eq. (IIV.9D . As* and v' 
contribute to the Higgs physics only as the combination of AeaV. Thus, four independent 
degrees of freedom are relevant in total. 

We fix the two parameters A$fc and AsaV in each of the following global fits for simplicity. 
Here, we consider the five possibilities, (A$fc, As^r/) = (1.0, v), ( 7 t,v), (0.1, v), (1.0, 5v) and 
(1.0, lOu) and search for their global minima on the remaining two variables ( m k ±±, sin a) at 
Erst. The positions of each corresponding minimum are (m k ±±, sin a) = (74.7 GeV, —0.0759), 
(109 GeV,-0.0774), (65.0 GeV, 0.102), (80.3 GeV,-0.0768) and (87.5 GeV,-0.0804) with 
Xmin — 5.08 (commonly in all the cases), respectively. The la and 2a boundaries are defined 
as the positions with y 2 = Xmm + A Xia, 2 a 5 where the values of Ay 2 (j2(r are calculated by 
the cumulative distribution function of the ^-distribution with two degrees of freedom as 
Axi a = 2.296, AyJo- = 6.180. The results are shown in Figs. HH5l Small disconnected 
regions are found possibly due to accidental cancellations when (A$fc, AeaV) = in 

Fig. [4] and (AA^V) = (1.0, lOu) in Fig. 0 When we consider the doubly charged scalars 
with m k ±± = 500 GeV, the range of a 


sin a) < 0.3 
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(IV.23) 















is roughly allowed within 2cr confidence levels, where actual upper bounds depend on contents 
of the two prefixed parameters. 


V. ISSUES ON N Rl DARK MATTER 


In this section, we have discussions on dark matter related issues of this model. After the 
breakdown of the new global 17(1) symmetry, an accidental Z 2 symmetry still remains, which 
ensures the existence of a dark matter candidate. Under the assignment of the 17(1) charge 
in Table [B Nr. (i = 1,2,3) and hf hold negative parities, where the lightest one among 
them gets to be absolutely stable. After considering the aptness discussed in Sec. 1IIIB1 that 
hf (and also hf) should be sufficiently heavy at around a few TeV to evade the constraints, 
e.g., via the processes with LFV, the lightest right-hand neutrino, namely N Rl , in our setup 
is usually stabilized by the symmetry and plays a significant role as dark matter. 

An important point of our dark sector is in the mechanism for generating Majorana mass 
terms of Nr that originates from the spontaneous global 17(1) breaking accompanying a 
nonzero VEV of the field So- Through the mixing between So and the Higgs doublet <3>, 
Nr 1 can communicate with the SM particles, where the two neutral scalars h and H work 
as mediators. Here, we expect that a suitable amount of dark matter relics are left in the 
Universe through resonant effects via the h or H pole, where the mass of N Rl is fixed around 
rrih/2 or m H /2 , respectively. 

The relic abundance and the spin-independent cross section for direct detection of the 
candidate N Rl are calculable by following a standard method. In the following part, we use 
the shorthand notation x f° r showing the DM Nr 1 (in the mass eigenstate as a Majorana 
fermion). Within the freeze-out approximation, the present-day relic density is evaluated 


as 


119] 


Q x h 2 ~ 


1.07 x 10 9 GeV -1 


(V.l) 


where h (~ 0.7), A7 P1 and g* express the present-day energy density of y, the present- 
day scaled Hubble parameter, the Planck mass and the number of the relativistic degrees of 


re latest value of £l x h 2 reported by the Planck experiment 
. The valuable Xf is defined by the dark matter mass m x 
and the temperature at the freeze-out Tf as m x /Tf. The efficiency of the annihilation after 


freedom, respectively. Note that t 
is 0.1196 ±0.0031 (68% C.L.) [120| 
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the freeze-out is described through the integral J: 


J&f) = 






(V.2) 


where (crw re i) stands for the thermally averaged annihilation cross section multiplied with 


the relativistic relative velocity u re i. In this work, we a c 


thermal average relativistically discussed in Refs. 


121 


opt 


the sophisticated way for taking 


122 ], where {av Te \) is estimated as 


ds- 


(av re i) = 


' 4m| 


^w xx iu (4 

16 xx V T 


' 4777-2 


ds x js — Arri^WxxKi ( —x 




miT 

\k 2 

2 


1 [K 2 (x)] 2 

X 

L V T J\ 





(V.3) 


where x is defined as m x /T (like xj) at the temperature T, Ad ,2 are the modified Bessel 
functions of the second kind of order 1 and 2, respectively, and W xx is a Lorentz invariant 
variable describing the cross section multiplied with the Lorentz invariant flux factor 4Ryu re i. 
W xx is formulated in the center of mass system with the integration over the solid angle as 


= E 


32vr 2 


1 - 


4m ^ f , _ 


7T\ I 2 


dSl\M(xx->ff)\ 


(V.4) 


where we sum over all possible two-body final states consisting of the same particle /. The 
detail of the amplitudes is found in Appendix O In the following numerical calculation, we 
simply use the fixed reasonable values Xf = 20 and g*( x f) — 100 throughout the calculations 
for brevity. 

Our e valu ation of the y-nucleon spin-independent cross section is based on the discussions 


in Refs. 


123, 


124], We consider the following effective Lagrangian for calculating the cross 


section at the leading order, 


» m cy 

C,s= E U™,XXm--£faXxG%G v “, (V.5) 

q=u,d,s 


where q, m q , a Sl and G° IJV represent the corresponding quark fields, the quark masses, the 
QCD coupling strength, and the field strength of the gluon, respectively. The coefficients 
fg (/ g ) determine the effective interactions between the quarks (gluon) and the DM x- The 
corresponding values in our model are 


u = 




cos a sin a. 


(V.6) 
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The spin-independent cross section with the target nucleus T is formulated by use of the 
spin-independent coupling of y with nucleon /tv ( N = p, n ) as 

<531 /^Ty I 'R'pfp T f^nfn\ j (V.7) 

71 A 


where / jlt x * s the reduced mass among the nucleus and the dark matter is defined as pr x = 
niT'm x /\rriT + m x ). tut shows the mass of the target nucleus T, which contains n p and n n 
numbers of protons (p) and neutrons (n), respectively. The effective y-nucleon (N) coupling 
/tv can be written down by the coefficients of the effective operators in £ eS (f q , f G ) and 
matrix elements (/wg, Ing) as 

/at/tuat = ^2 fqfwq + g/c/vc, (v.8) 

q=u,d,s 

with the concrete forms of the matrix elements, 


(N | m q qq\N) = f Nq m N: 1 - ^ fif q = Ing, 

q=u,d,s 


(V.9) 


where the value of /tvg is calculable through the latter relation by the use of fwg- Like in 


Ref. 124], we adopt the following default values in the program micrOMEGAs 125]: f nu = 
0.0110, f nd = 0.0273, f ns = 0.0447. 

The latest bound on the spin-independent scattering process was reported by the LUX 
experiment as an upper limit on the spin-independent (elastic) DM-nucleon cross section, 
which is approximately 10 -45 cm 2 (when m x ~ 10 2 GeV) with the 90% confidence level 126 ], 
Here, we choose the neutron as the nucleon for putting a bound on our parameter space 
via the LUX result. Now, the (spin-independent) y-neutron cross section is calculated with 
ease as 

2 


a, 


N=n 

SI 


/Gy 
71 A 


m r 


^ fnq T fnG j fq 
q=u,d,s / 


z z s~i Z 
= -^n X m n C 


2 2 ( m X COS a S111 a 


7T 


vv' 


-2 + 

mi 


m 


H 


(V.10) 


with 


c = l + l X /»«- 0 ' 287 . 

q=u,d,s 

where similar calculations were done, e.g., in Refs. 


(V.ll) 


127 


128|. 
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FIG. 6: Realized present-day relic densities of the dark matter candidate x are shown as functions of 
m x . The other effective parameters are fixed as sin a = 0.4, v' = 800 GeV, mu = 250 GeV, m k ±± = 
500 GeV (left panel) and sin a = 0.5, v' = 5TeV, mu = 250 GeV, m k ±± = 500 GeV (right panel) 
for demonstration purposes, respectively. The area inside the two horizontal dashed lines (where 
the splitting is al most invisible) suggests the 2<r consistent region with the Planck data [0.1196 ± 


0.0031 (68% C.L.) 


1201 ]]. In the current setup, the relic density sharply drops around the two 


resonant regions around m^/2 ~ 62.5 GeV and mjj /2 — 125 GeV, where in the latter case, shown 
in the right panel, dropping is not sufficient for obtaining a proper amount of dark matter even 
around m#/2 at present. We neglect the three- and four-body final states via virtual W and Z 
boson decays, which gives sizable modifications near the thresholds for producing gauge boson 


pairs 


129), 


1301 ]. since our interest is only around rn k /2 and mn /2 where this correction is expected 


to be insignificant. 

In the following calculation, we treat the variables a, m #, and v' independently, which 
determines the coefficients A®s, A®, and As through the relations in Eqs. (III.61) (III.71) as 


A<j>s — 


A® — 


As = 


sin a cos a m 


H 


m 


vv 

cos 2 a m 2 + sin 2 a m 2 H 
sin 2 a m 2 + cos 2 a m? H 


2v' 2 


(V.12) 


We set two quartic couplings A®*, and As k as 1.0 and 0.1. In our scenario, v' tends to 
be 0(1) TeV leading to suppression of the thermally averaged cross section describing pair 
annihilation processes; see Appendix O for detail forms. It would require a mechanism 
for enhancing the cross section. Hence, in this manuscript, we focus on the two resonant 
regions of m x around m^/2 or mu/ 2. Our requirement for the amount of the relics is that 
it should be within the ±2<r range of the latest value reported by the Planck experiment 
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FIG. 7: The matrix plots indicate suitable choices of the two parameters v' and sin a when we 
consider the situation of mu = 250 GeV, m k ±± = 250 GeV. We obtain the observed relic density 
in the red (blue) region in the left (right) panel when m x is around m/,/2 ~ 62.5 GeV [mu I 2 — 
125 GeV), respectively. The green region is excluded via excess of the invisible decay channel of the 
observed Higgs boson. No excluded region is found by the direct detection in the shown parameter 
range. 


v' (TeV) 


v' (TeV) 




FIG. 8: The matrix plots for showing the region where a proper amount of the relic density is 
realized when mu = 250 GeV, m k ±± = 500 GeV. Conventions are the same with ones in Fig. (7J 
No excluded region is found by the direct detection in the shown parameter range. 


[0.1196 ± 0.0031 (68%C.L.) 120]]. In the following matrix plots, the values of v' and sin a 
at the center of a square represent the inputs in the whole region shown by the square. 

An important point is that all the effective diagrams of this pair annihilation of \ are Higgs 
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FIG. 9: Similar plots as Figs. [7] and [8] when m# = 500 GeV, m k ±± = 250 GeV (left) and rriu = 
500 GeV, m k ±± = 500GeV (right), respectively. Here, no solution is found around m x — m#/2 ~ 
250 GeV, while one around m x — rrih/2 ~ 62.5 GeV is still there. The green region is excluded by 
the invisible decay of the 125 GeV Higgs. No excluded region is found by the direct detection in 
the shown parameter range. 


mediated. 8 Like the couplings between the observed fermions and the 125 GeV Higgs boson 
in the SM, the DM-portal scalar couplings are proportional to the factor m x /v'. In Fig. [6] 
we present realized present-day relic densities of the dark matter candidate x as functions 
of m x , where the other effective parameters are fixed as sin a = 0.4, v' = 800 GeV, mu = 
250 GeV, m k ±± = 500 GeV (left panel) and sin a = 0.5, v' — 5 TeV, mu = 250 GeV, m k ±± = 
500 GeV (right panel), respectively. Here, we can see that a suppression due to a huge v' in 
the latter case and the point around the heavier resonance of H is not suitable for a correct 
amount of DM today. Around the lighter one of h, less s-channel propagator suppression 
could be expected and thus the solution around this pole still remains. 

Another critical variable in the relic calculation is the mixing angle a among the two 
CP even scalars h and H. If sin a is close to zero, contributions from the mediator h are 
diminished significantly due to almost zeroness of the couplings between x and the SM 
particles in the final state. Thereby, a suitable amount of the mixing is requested for the 
solution near m^/2. In the following analysis, we focus on non-negative values in sin a 


Note that the processes 
important in similar situations 


—> itt ■ (i,j = 1,2,3) with exchanging hf in t-channel, which are 
13l|, are ineffective and neglected in the present scenario because the 


lower bound on m h ± is around a few TeV and decoupled as shown in subsection IIII Bl 


34 



































































































since the sign of sin a is insignificant in the resonant regions. 9 Also, we ignore the range of 
| sin a\ > 0.5 since no allowed region was found in the results of the global analyses shown 
in Figs. 0H5J 

As actual examples, we investigate the four possibilities in (rriH,'n r ik±±) of 
(250 GeV, 250 GeV) in Fig. [71 (250 GeV, 500 GeV) in Fig. [HI and (500 GeV, 250 GeV) and 
(500 GeV, 500 GeV) in Fig. [9l respectively. The regions covered with the color green are 
excluded, where the branching ratio of the observed Higgs boson to invisible pairs exceeds 
the 95% confidence level upper bound 0.29 reported by the ATLAS experiment |132j |. 10 The 
value of v' and a almost completely determines the profiles since the primary final state 
with invisibleness is a pair of the NG bosons G. We mention that no exclusion from the 
direct detection process is found in the shown parameter ranges since the elastic scattering 
cross section is suppressed by the large VEV v'. We observe small differences originating 
from the value of rrik±± since in the DM mass range of m x < rrik±±, fc ±=t can contribute to 
the relic density only indirectly through the total widths and T#. 11 

A suggestive aspect via the results is that only around the m^/2 is possible when rrin gets 
large as 500 GeV because the suppression via the heavy resonant particle could eliminate 
the solution near m#/2. On the other hand, as discussed in Sec. I1V1 the 7 and 8TeV LHC 
results restrict the scalar mixing angle a , where we quote a typical bound we obtained, 

| sin a| < 0.3. For making the DM solution around rxih )2 realizable, we should select the 
absolute value of the mixing angle | sin a\ around or more than 0.3 as shown in Figs. [7]-[9j Let 
us remind the reader that the sign of sin a does not affect the result (calculated in sin a > 0) 
significantly. In addition, here, the value of v' is also constrained as around or less than 
9TeV in the case of m x ~ rrih/2. In conclusion, the present scenario is still viable, but a 


9 We mention details of this point. The sign of sin a affects results only through the amplitude of XX hh 
and the partial decay width of h/H —> 77 and h/H —>■ Z 7 in the s-channel propagators. In the DM 
masses which we consider, the amplitude gives only a subleading contribution. The effect coming from 
the widths is also negligible since the branching ratios of h/H —>■ 77 and h/H —> Z 7 are pretty small. 

10 In this estimation, the Standard Model production cross section was assumed in mu = 125 GeV. Here, 
we do not care about this point seriously and simply utilize this result for putting a bound on the present 
scenario. 

11 We simply ignore the process XX k ++ k in the case near the kinetic boundary, mu = 500 GeV and 
m^±± = 250 GeV, where the heavier resonant region corresponds to m x ~ 250 GeV. Here, we do not 
consider the lighter k ±:t with 100 ~ 250 GeV mass, the possibility of which was examined in the global 
fits in section m 
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part of the parameter space generating a suitable amount of relics comes to be on the edge. 
Note that our calculation would be applicable for a rather general setup of a fermion dark 
matter communicating the SM sector through two scalars with renormalizable interactions 
and 0(1) TeV v' after the decoupling of fc ±=t . 


VI. SUMMARY AND CONCLUSION 


In this paper, we proposed a three loop radiative neutrino mass scenario with an isolated 
doubly charged singlet scalar k ±zL without couplings to the charged leptons at the leading 
order for solving an undesirable situation in our previous study [l|. The previous scenario 
has the same particle contents of charged singlet scalars with the present model, two singly 
charged ones ( hf , h^) and one doubly charged one (fc ±=t ), but assignments of additional 
global symmetries are different, which makes a difference in interactions associated with the 
charged bosons. 

In the previous case without Majorana neutrinos, k ±zL should attach the charged leptons 
at the leading order to generate a neutrino mass matrix at the three loop level. But a side 
effect due to this type of couplings with LFV could be serious from a phenomenological point 
of view. Because of resultant tree-level contributions to LFV processes, related couplings 
are severely constrained, where realized neutrino masses get to be smaller than the observed 
values. Increasing scalar trilinear couplings with charged scalars appearing in the numerator 
of elements of the loop-induced neutrino mass matrix seems to be an option. However, the 
largeness in the trilinear couplings tends to destabilize the vacuum rapidly. Making all of the 
charged scalars very heavy is a solution since the scalar trilinear couplings do not contribute 
to LFV at the leading order and the ratio of the trilinear couplings divided by the masses of 
the charged scalars can remain “sizable” for the observed neutrino masses. Here, however, a 
typical mass scale of the charged bosons should be around 10 TeV, where we cannot expect 
a signal suggesting the existence of the particles at the LHC. 

Our modified setup in the assignment of a global U( 1) symmetry brings us a better 
situation. Introducing Majorana neutrinos allows /c ±=t to be away from the charged leptons. 
Now the LFV is relatively relaxed, and at least a part of the couplings contributing to 
the neutrino mass can take larger values. Thus, scalar trilinear couplings need not to be 
drastically large, and, as a result, vacuum destabilization is alleviated. We found that a 
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few hundred GeV /c ±± is realizable, which can be tested at the ongoing LHC experiment, 
even though the two singly charged scalars are still very heavy at around a few TeV. It is 
noteworthy that the associated NG boson G does not generate dangerous interactions that 
severely restrict the decay constant v'. The constant v' can takes a value around the TeV 
scale with no harm. 

Another constraint on the system comes from the mixing of two CP even neutral scalars, 
which are mixed states between the component of the SU(2)l doublet $ in the SM and the 
corresponding part of a complex SU(2)l singlet scalar S 0 for the spontaneous symmetry 
breaking of the global U(l). After the breakdown, Majorana fermions obtain masses via the 
Higgs mechanism and a residual discrete parity ensures the stability of an accidental DM 
candidate of the lightest Majorana neutrino N^. This mixing angle plays a very important 
role in the DM-related phenomena in this model since the candidate can couple to the SM 
sector only through the two mass eigenstates of the CP even scalars h (= 125 GeV) and H . 
Note that the two singly charged scalars should be decoupled in our scenario and diagrams 
associated with them are negligible. 

The latest LHC result of the Higgs search by the ATLAS and the CMS experiments 
puts a stringent bound on the mixing angle a, | sina| < 0.3 as a roughly estimated bound 
through a global analysis. A typical scale of the VEV of S 0 would be around a few TeV scale, 
where the DM pair annihilation cross section is significantly suppressed since the Yukawa 
couplings between Nr ± and So are inversely proportional to the VEV. When the DM mass 
is around the two resonant regions, m/,,/2 and m#/2, the cross section is enhanced very 
much. But we found that when Mh is heavy at around 500 GeV, the possibility near m#/2 
is closed due to the propagator suppression of the massive mediator. On the other hand, to 
activate the solution around m/,/2, the mixing angle a should be large to a certain degree to 
ensure a sizable interaction between the DM and SM particles. Typically, a required value 
in this direction is |sina| > 0.3, which would not have large overlaps with the preferred 
area obtained via the global analysis on the Higgs search results. 

In spite of the insignificance of the bounds via the invisible 125 GeV Higgs decay and 
the DM direct detection because of the coupling suppression via the huge VEV of E 0 , 
experimental data in Higgs physics put sizable constraints on our scenario. In the near 
future, updated Higgs results could declare validity of our scenario much more precisely. 
Also, searching for a suitable way to discriminate the present scenario from other models 
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with charged scalars will be an important task at the LHC and other future colliders. 
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Appendix A: Analytic forms of loop functions for LFV 


as 


In this appendix, we show analytic forms of the loop functions for processes with LFV. 
The functions for l 1 7 in Eq. (1II.34|) are obtained 
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Note that when we take the limit M^ a —> 0, I[ a and I' 2 a are reduced to /i iQ and l 2 , a (up to 
the difference of m h ± and m h ±), respectively, 
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The functions J\. a b and J 2 , a b for describing the £ —> 3£ processes in Eq. (111.40j) take the 
following forms 
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where the following specific cases are also obtained: 
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Appendix B: Partial widths of CP even scalars 


In this appendix, we represent a part of partial decay widths with nontrivial forms of the 
two CP even scalars h and H with our notation for loop functions. 
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with Chkk in Eq. (IIV.OI) and the parameters 
, 4m 2 


T; = 


m 


(i = t, W, k), C Hkk = sin aX$ k + cos aA Sfc — . 


H 


(B.12) 


G F and a s are the Fermi constant and the QCD coupling strength, respectively. Note that 
in the calculation Vh^zj, the replacement 7/ —* r- in the factor A^ is required. One can 
refer to Sec. IIVI for the detail of loop functions. 
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Appendix C: Matrix elements for N Rl dark matter annihilations 


In this appendix, we summarize the averaged matrix elements squared for relic density 
calculation of the dark matter candidate N Rl through pair annihilations. Note that no 
symmetric factor is included originating from identical final states in the following forms, 
which should be considered in the integration over the solid angle for estimating the thermal 
average of the cross section (cru re i). 
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\M(N Rl N Rl hh)\ 2 = 


(m N i \ 2 
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s-m 2 h + im h T h 

s — m 2 H + imn^H 
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[(pi -P 2 ) -Ml J , 


(C.4) 


(C.5) 


where Nc f shows the color factor 3 (1) when / is a quark (lepton), pi,p 2 and k\,k 2 are 
initial- and final-state momenta, respectively, s is defined as s = (pi + P 2) 2 = {k\ + k 2 ) 2 . 
Here, we use the shorthand notations of c a = cos a and s a = sin a. The forms of effective 
couplings describing scalar self interactions are 


C hhh = 3A$ S (v(s 2 a c a ) - v\c 2 a s a )) +6A $vc 3 a - 6A xv's 3 a , (C.6) 

C H hh = A<j, s (w (4 _ 2 4 s a) + v' (4 - 2 c a s 2 a )) + 6A q,vc 2 a s a + 6A ^v's 2 a c a . (C.7) 

We mention that the processes N Rl N Rl into HH, hH , k ++ k , and are inef¬ 

fective in our discussion when the mass of the DM is around m^/2 or m R / 2 considering the 
bounds discussed in Sec. mm 
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Note that the processes —>■ £fL (■ i,j = 1.2, 3) with the exchange of h f in the t 


channel, which are important in similar situations 


7, 


131 ], are ineffective and neglected 


m 


the present scenario because the lower bound on m h ± is around a few TeV and decoupled 
as shown in Sec. IIIIB1 We neglect the three- and four-body final states via virtual W and 
Z boson de cays, which gives sizable modifications near the thresholds for producing gauge 
boson pairs 129, 130], since our interest is only around mh /2 and m#/2 where this correction 
is expected to be ineffective. 
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